Deterministic coupling of a single nitrogen vacancy center to a photonic crystal cavity 
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We describe and experimentally demonstrate a technique for deterministic coupling between a 

photonic crystal (PC) nanocavity and single emitters. The technique is based on in-situ scanning 

of a PC cavity over a sample and allows the positioning of the cavity over a desired emitter with 

nanoscale resolution. The power of the technique, which we term a Scanning Cavity Microscope 

^^ (SCM), is demonstrated by coupling the PC nanocavity to a single nitrogen vacancy (NV) center 

JT^ in diamond, an emitter system that provides optically accessible electron and nuclear spin qubits. 

(N 

^.^ Optical resonators enable large amplification of small optical signals, resulting in a range of spectroscopic and 

C^ sensing applications, and have allowed for detection of single atoms [1, molecules [2 , and quantum dots [3l[4]. In 

^) addition, they enable a controllable coupling between optical emitters and the cavity vacuum field that is critical for 

^^ efficient light sources [5]-[7] and for the realization of memory nodes in quantum networks [8 and quantum repeaters 

CN [9]. This coupling strength scales with the cavity mode volume Vm as 1/v^T^, and consequently, nanoscale photonic 

^"^ crystal (PC) cavities have been explored extensively in solid-state cavity QED applications. While much progress 

^ has been achieved in coupling quantum dots to PC cavities made from the host material [3,|T0l[TT], extending these 

r"| techniques to fully deterministic coupling and to other material systems has been difficult. Specifically, there has 

Qh been much recent interest in coupling PC resonators to NV centers [T2]-[15], a promising single photon emitter with 

4L excellent electronic and nuclear spin memory pIBHIS] . though experimental demonstrations have remained a challenge. 

rt We demonstrate a technique for deterministic positioning of micron-scale PC slabs that support high quality factor 

(Q) cavity modes with nanometer-scale features. When such a cavity is scanned over the sample, it can be used 

for deterministically coupling to optically active systems with sub-wavelength resolution via the evanescent field. By 

appropriate design of PC cavities and waveguides, these systems combine sub- wavelength resolution, high throughput, 

and cavity-enhanced sensitivity. In particular, they can be deterministically interfaced with isolated optical emitters. 

In our experiments, the PC consists of a triangular lattice of air holes in a gallium phosphide (GaP) membrane, 

creating an optical bandgap that confines light in the slab to a cavity region. The bandgap along the TX crystal 

direction is shown in the dispersion diagram in Figure Ilia). Confinement in the vertical direction occurs through 

^N^ total internal reflection (TIR) for modes with frequencies below the air light-line indicated in Figure Ilia). A row 

(SJ of missing holes supports band modes which form bound cavity states when terminated on two sides. We employ a 

^jJ three-hole defect cavity [19] whose geometry is optimized for use on a Poly(methyl methacrylate) (PMMA) substrate 

^^ with a refractive index of n^ ^ 1.5. The TIR-confined region in k-space is smaller on top of the PMMA, as sketched 

^^ in Figure ma), but simulations indicate that the Q value can still be above 13 • 10^. The cavity has a mode volume 

^— H Vm = 0.74(A/nGap)^, where ncaP = 3.4 is the refractive index of GaP at A = 670 nm. The fundamental mode of the 

K^ PC cavity is depicted by its energy density in Figure [ijc). The cross section in Figure [ijb) shows the evanescent tail 

. ^^ of the mode that couples to emitters. 

^\^ We fabricate GaP PC nanocavities by a combination of electron beam lithography and dry etching [20] of a 108 

J-H nm membrane of GaP on top of a 940 nm-thick sacrificial layer of a Alo.85Gao.15P. A wet etch removes the sacrificial 

layer, leaving free-standing photonic crystal membranes. The scanning electron micrograph (SEM) of a resulting PC 
nanocavity is shown in Figurellld). Reflectivity measurements of freestanding cavities indicate that quality factors (Q) 
of these cavities can exceed 6 • 10^, the maximum value that can be measured with the resolution of our spectrometer 
(Figure [ijf)). However, in the remainder of this paper we will study cavities with typical Q values below 1000, since 
these were more reliably fabricated in large numbers, permitting systematic studies. To transfer cavities, we press the 
GaP chip against a flexible polymer layer of Polydimethylsiloxane (PDMS), which separates the PC membranes from 
the chip while preserving their arrangement. The adhesion between the membranes and the PDMS is weak enough 
so that the GaP structures can be stamped onto the sample that is to be imaged, as shown in Figure [ije). In our 
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FIG. 1: (a) Dispersion for the photonic crystal slab in air, along the waveguide direction kx] the inset shows the crystal (orange) and 
inverse crystal (blue) directions. The lattice has a periodicity of a = 176 nm, hole radius of ~ 53 nm, and slab height of ~ 110 nm. (b) 
Energy density for fundamental mode in cross-section and (c) in plane, (d) SEM. (e) The photonic crystals are transferred from the GaP 
chip onto a substrate via a polymer stamp, (f) Broad-band reflectivity measurement of a cavity resonance with Q ^ 6 • 10^. (g) The PC 
slab is positioned relative to a target nanocrystal in the polymer fllm. 



demonstration, the sample consists of ^ 30 nm diamond nanocrystals that are dispersed on a glass slide covered by 
a 100 nm thick layer of PMMA, for which the transfer process succeeds with ~ 80% probability for each membrane. 
The sample is mounted in a scanning confocal microscope with an oil immersion lens. A tungsten tip with radius 
< 0.5 /im is used to scan and position a PC nanocavity with nanometer resolution (Figure [ijg)). 

Figure [2|a) shows a room-temperature photoluminescence (PL) image of a typical sample (red signal) obtained by 
scanning a green excitation laser across the surface with galvanometric mirrors. The PL spots in the image correspond 
to single or clusters of NV centers in diamond nanocrystals. The PC nanocavity can be located with the excitation 
laser reflection (green in Figure [2|a)) as well as weak fluorescence originating from impurities in GaP and PMMA 
(green curve in Figure |2|c)). The PC fluorescence clearly shows two resonances at Ai = 667.3 nm and A2 = 643.0 
with quality factors Qi = 550 and Q2 = 610. These quality factors are rather low because of variabilities in the 
fabrication; we saw no degradation due to the transfer onto the PMMA substrate. By collecting spectra at different 
points within the cavity, we can identify peak 1 as the fundamental mode depicted in Figure [ijc) and peak 2 as two 
nearly degenerate, oppositely polarized higher-order modes of the cavity. The technique and results of the spatially 
resolved cavity spectrum are discussed in Appendix [Cj 

We achieve deterministic coupling between the NV and the nanocavity by flrst selecting a 'target' NV (indicated 
in Figure [2| a)). This center exhibits a broad spectrum /q (Figure [2|c)), which is characteristic of NV centers and 
results from a broad phonon sideband extending from ~ 640 to 800 nm. Importantly, the emission exhibits a strongly 
anti-bunched autocorrelation (inset), indicating that it results from a single emitter. To couple this emitter to the 
cavity, we position the PC membrane over the target NV using the tungsten tip. As shown in Figure [2|d), the PL 
spectrum Ic changes dramatically and shows strong peaks on resonance with the cavity modes. The intensities of 
these peaks are far higher than the cavity background, Icb- Moreover, the autocorrelation of the coupled NV-cavity 
system is again strongly anti-bunched (Figure |2Vd)), indicating that it is driven by emission of the NV center. We also 
verifled the electronic triplet state of the coupled NV by electron spin resonance measurements, as shown in Appendix 



The SE rate of an NV center is also modifled by the presence of the PC slab. Speciflcally, Figures [2|e) and (f) 
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FIG. 2: The photonic crystal is moved from an initial uncoupled position (a) into alignment with the target NV center (b). The pump 
laser reflectivity is shown in green and the photoluminescence in red; pump laser power is 500 /iW, focused to ~ 0.2 /xm. (c) PL spectrum of 
the uncoupled NV (/q) and uncoupled cavity background (Icb)- A photon correlation measurement (see Appendix IbI shows that the NV 
emission is strongly antibunched (inset); this feature is surrounded by photon bunching due to shelving in a metastable state of the NV 
emitter [22 . (d) PL spectrum Ic of the coupled NV-cavity system, again strongly antibunched (inset). A flt to theory (Eq[l| gives the SE 
rate into the cavity normalized by the background emission rate, f^(X2) = 5.3, /^(Ai) = 0.7. (e) Time-resolved emission for the uncoupled 
NV, far removed from the PC membrane, and (f ) the coupled NV. The 6 ps excitation pulse was generated by a frequency-doubled 1064 
nm laser at 20 MHz repetition. 



show that the hfetimes of the uncoupled and coupled NV centers are ro,c = 16.4 ± 1.1, 12.7 ± 1.5 ns, respectively. The 
PL spectra on and off the PC coupled with lifetime measurements allow the determination of the spectrally resolved 
SE rate enhancement, F(A), of the coupled emitter via the relation F{X) = Ic{\)tq/Iq{\)tc (see Appendix [A|) : the 
analysis of the data in Figure |2] yields i^(Ai) = 2.2 and F{\2) ~ 7.0 (the full curve F{\) is plotted in the Appendix). 

We next demonstrate the spatial resolution of our method. By monitoring the fluorescence spectrum while scanning 
the cavity over the NV, we can map out the near- field emitter-cavity coupling. This is demonstrated in Fig 3(a-e), 
where we scan the cavity along its longitudinal (x-axis) over the sample in 3.4 nm steps. Figure 3(f) presents a 
series of PL spectra acquired as the cavity moves over the emitter, and reveals an intensity oscillation with a period 
corresponding to one PC lattice spacing, a ~ 180 nm. This oscillation corresponds to the spatially dependent SE 
modification, which is directly proportional to the cavity's electric field intensity. 

To analyze our observations, we note that the fluorescence of the coupled NV-cavity system is given by the emission 
directly from the NV, the emission through the cavity, and interference between the two: 



5d(w,r-) = Cnv + Ce„„r(r-)|L(w)|2 + 2C„tSR[e^^*\/7^L(w)], 



(1) 



where Cnvi Ccavi and Cint determine the relative contributions of the NV, the cavity, and their interference, re- 
spectively, which depend on the collection geometry and coupling to the collection fiber. This relation is derived in 
Appendix If] L{uj) = 1/(1 + i{uj — uJc)/k,) gives the Lorentzian line shape of the cavity resonance at ujc with linewidth 
hi = uJc/^Q^ and Acj) accounts for the phase difference at the collection point between the direct NV emission and the 
emission through the cavity. The factor /^(cj, r) is the SE rate enhancement of transitions in the phonon sideband of 
the NV with respect to the background emission rate into non-cavity modes. 

The coefficients Cnv^ Ccavi and Cint can be estimated from our experimental data as follows. Because of the 
high numerical aperture of our objective, nearly half of the emission from the cavity and the NV is collected: this 
observation suggests Cnv ^ Ccav When the signal is collected through a single-mode fiber, the interference term 
represented by Cint becomes important and results in Fano-like features in the spectrum (see Appendix |d| [22] . 
However, we find that the interference term vanishes when a multi-mode fiber is used, and we can set Cint =0- A fit 
of EqjTlto the spectrum in Figure ^d) then yields /^(A = 643 nm, r) = 5.3, /'^(667nm, r) = 0.7. 
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FIG. 3: Scanning of the PC nanocavity probe in small steps, shown in snap-shots (a-e). (f) Photoluminescence scans for (Ax) = 3.4 nm 
average step sizes, (g) Fitted cavity SE rate enhancements /^(Ai, r) for mode 1 showing a FWHM resolution of A ~ 80 nm. (h) Expected 
SE enhancement factor /^(Ai, r) and the estimated trajectory of the NV at Az = 98 ± 5 nm, Ay = 70 ± 5 nm, and jl in the plane at 20° 
to the X-axis. The indicated track matches the observed SE enhancement in (g). 



Since the signal in Figurejsj^f) is proportional to Sd{oj^ r), we can now use EqjTJto compare the measured cavity signal 
to theory. Figure IsFg) plots the fitted values of /^(cji,r) for the fundamental cavity mode frequency uji = 27rc/Ai, 
as shown in the red crosses. By comparing the experimental f^{uji^r) values to predictions for the cavity mode, we 
find a match between experiment and theory for an NV dipole /i that is z = 98 ± 5 nm from the PC surface, as 
expected from the PMMA thickness, and at an angle of 20° to the x-axis, obtained from the best fit to the data. 
For these conditions, the predicted value of SE rate modification corresponds to the track graphed in Figure |3|h), 
in good agreement with experimental observations. A small discrepancy in the fit at Ax ~ 190 nm results primarily 
from positional slip of the PC cavity that can build up during the scan, a problem which could be improved by rigidly 
attaching the membrane to a stiffer scanning tip. 

The high spatial resolution and frequency- selective modification of spontaneous emission opens new possibilities for 
efficient interfacing of promising solid state qubits via optical fields. For instance, while the NV center is a promising 
system for quantum information processing, only the emission occurring into the zero phonon line (ZPL) is suitable for 
coherent optical manipulation. The frequency-selective emission enhancement demonstrated here potentially allows us 
to direct most of the emission of the selected NV centers into the ZPL. Furthermore, the hybrid approach is compatible 
with narrow linewidth NV emitters in bulk diamond at low temperature. This opens the door for applications ranging 
from quantum repeaters to single photon nonlinear optics. Moreover, although we have focused here on NV centers, 
our scanning technique provides a 'cavity QED interface' that can be of use to a broad range of solid state qubits. 

Furthermore, the PC scanning technique can serve as a new imaging approach with sub-wavelength resolution and 
high throughput, which we term a Scanning Cavity Microscope (SCM). Unlike other near-field probes that compromise 
the signal intensity to achieve high spatial resolution, SCM enables large count rates: in the demonstration shown 
here, we record up to ~ 1 • 10^ photons/s from a single NV, exceeding collection rates reported to date with far field 
optics. This can be further improved by efficiently out-coupling through cavity-coupled waveguides. In addition, the 
spatial resolution of the SCM is determined by the feature size of the confined field, which is A ~ 80 nm for this 
cavity. This in-plane resolution may be improved substantially using cavity modes with small feature sizes, as in slot- 



waveguide cavities [23]. These qualities make the SCM a promising tool for label- free single molecule studies [21124] 
or high-resolution studies of local index variations in thin films [25j [26] . We note that a scan acquired by moving the 
photonic crystal membrane across the surface may need to be de-convolved by the response function of the cavity 
mode, as described in Appendix [Gj Beyond high resolution and throughput, the SCM adds the capability to modify 
the spontaneous emission rate to near-field microscopy. This opens new possibilities for direct investigations of decay 
channels of optical emitters, such as light emitting diodes or fluorophores; for instance, by monitoring the emission 
intensity while effecting a known change in the radiative emission rate, the relative nonradiative recombation rate 
may be inferred, allowing a direct estimate of the radiative quantum efficiency of the material. 

Appendix A: Wavelength- resolved excited state recombination rates 

The excited state of the isolated NV~ decays through a set of Nx radiative channels which we distinguish by their 
wavelength Xj and describe by the recombination rate 7^, and nonradiative recombination described by a term Jnr- 
We consider the NV at some position fi from the cavity, and assume that it is pumped into the excited state at some 
rate pi. Then the emission spectrum at wavelength j is given by 

lij = cpiTjij S^\^^ — (Al) 

= cp,r]^J^ (A2) 

where r]ij is the collection efficiency, F^ = ^^ 7^/^ + jnr is the excited state decay rate, and c is a proportionality 
constant. The term I3ij = ^ij/Ti can be identified as the fraction of the total recombination rate that occurs through 
wavelength j. The decay rate can be obtained from a lifetime measurement, F^ = l/r^. The ratio of the spectra at 
position i and the 'bulk position' (away from the PC slab) gives: 

J-ij _ PiVijlij^O ( A'^) 



If the pump rate and collection efficiency are equal, we have 

J-ij Tij-'- 

To 



(A4) 



= F,,-^, (A5) 

J- i 

where the Purcell factor Fij = Jij/joj = lijTi/IojTo. For a position i = 1 where the NV and cavity are aligned on 
top of each other, Fij = F{Xj) is plotted in Figure Is] 

Appendix B: Photon correlation measurements 

The photon correlation measurement is performed with a Hanbury Brown / Twiss (HBT) setup. The PL is directed 
to a beam splitter, after which the two arms are coupled through single mode fibers (SMF) to single photon detector 
modules (Perkin Elmer). A time-correlated counting module is used to assemble the histogram of clicks on the two 
detectors. Figure |4] shows the spectra collected through the SMF for the coupled and uncoupled positions of the 
cavity. 

Appendix C: Characterization of cavity modes 

To further characterize the cavity resonances experimentally, we measured the polarization dependence of the 
emission collected from two different locations while pumping the NV near the center of the cavity. These positions 
are indicated as xi and X2 in Figure |6ja), together with the simulated electric field energy density for three relevant 
cavity modes. Collecting from the pump location X2, we observe the mode at 667 nm to be strongly polarized, while 
the mode at 643 nm exhibits only weak polarization dependence (Figure p[b)), for a total of three modes. From xi, 
we only observe polarized modes at Ai and A2. From these observations, we conclude that the peak at A2 corresponds 
to the two modes shown in Figure|6ja.l,a.2) with nearly degenerate frequency (within the relatively broad linewidth), 
opposite polarization, and different spatial extent. The peak at Ai indeed corresponds to a single mode, which we 
identify as the lowest- frequency confined state of the L3 cavity. 
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FIG. 4: Spectra collected into single mode fiber, (a) PL from uncoupled NV and cavity, (b) PL from cavity-coupled NV centre, (c) 
Scanning the cavity over the NV centre produces interference effects between NV and cavity emission. Fits to Eq.l with C2 = 0.6. 
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FIG. 5: The calculated spontaneous emission rate modification F{X). 
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FIG. 6: Polarization dependence of cavity emission collected from spatially distinct points. We position the cavity so that it slightly 
off-center from the NV emitter. Exciting the NV and collecting from the point of excitation (a) we observe the fluorescence at 667nm to 
be strongly polarized while the fluorescence at 643 nm is only slightly polarized. However, by collecting from the end of the cavity in (b-c), 
we can pick out the resonance at 643nm which is polarized along the same axis as the 667nm mode. 

Appendix D: Single mode fiber spectra 

When the multi-mode fiber was replaced with a single mode fiber, the spectra in Fig. 2(c,d) in the main text show 
the same pattern, but at a reduced count rate, as displayed in Figure [4J 

The simultaneous collection from the cavity and the emitter results in interference at the detector. As the cavity is 
scanned over the NV emitter, we observe Fano-like features at the cavity resonances shown in Figure [4Fb), which are 
fit well by Eq.l with C2 = 0.6. Such interface effects have been described before in microdisks [22] and microspheres [27] 
coupled to single NV emitters. In microscopy, the intereference signature may be used to better distinguish a single 
emitter from the background. 



Appendix E: Electron Spin Resonance and Rabi Oscillations 



The NV center has two unpaired electrons in a triplet configuration with a zero-field splitting of A = 2.87GHz 
between m^ = and m^ = ±1 sublevels. The excited state decay rates allow spin-sensitive detection since the 
fiuorescence is reduced for the mg = ±1 states [28 . In the spin measurements described here, the mg = ±1 sublevels 
were additionally split by an external magnetic field produced by a permanent magnet. Figure iTFa) shows the 
microwave transitions from rris = to vfis = ±1, which are evident in the intensity of the cavity-coupled NV as a 
function of a microwave field of intensity v^ while Figure JTlb) shows driven spin oscillations when the NV spin is 
optically initialized in m^ = 0, excited with a microwave pulse for a duration t at a frequency v = 2.77 GHz, and 
optically detected using a pulse sequence described in elsewhere [29] . 
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FIG. 7: (a) Spin transition observed in the NV fluorescence under excitation with a microwave fleld at frequency u. (b) Driven spin 
oscihations for u = 2.77 GHz. 



Appendix F: Spontaneous emission modification 

The Hamiltonian describing a NV (lowering operator a = \g) {e\) coupled to a cavity mode (with annihilation 
operator a, without considering environment interactions, is given by 



H = —a^a - -^{cTz) + ig{cra^ - aa^), 

where A is NV-cavity detuning and the reference energy is the mean of NV and cavity energy. 
The Master equation is given by 



dp 



-i[H,p] + -{2apa^ - a^ap - pa^a) + -{2a pa^ - a^ap - pa^a) + —-{(JzpcJz - p) 



(Fl) 



(F2) 



at, a]. 



where 7, n and 7^^ account for NV population decay, cavity population decay and NV pure dephasing; and <j. 

We co nsider the subspace of or 1 excitation in the NV \g) , |e), and or 1 excitations in the cavity |0) , |1) . From 
Eq. F2 we obtain the Maxwell Bloch equations |3Q] 



I da 



da\ , . A K.. , , / . 

A 7 



(F3) 

(F4) 
(F5) 



These two coupled equations are solved with the initial condition a(0) = 0,<j(0) = 1; {(izCi) = —a in the subspace 
considered here. We define Ci = {—i^ ~ f ) ^^^ ^2 = {^^ ~ 2 ~ ^d) and A_ = Ci + C2 — \/(ci — 02)^ — 4^^ , 
A+ = ci + C2 + \/{ci — ^2)^ — 4^^- The time-dependent solution is 

e 



a{t) = 



j^gA+t _gA_tj^ 



and 



a{t) 



\/(ci - C2)2 - 452 

(ci - C2)(e^-* - e^+*) + \/(c7^^c^F^^V(e^-* + e^+*) 



2V(ci - C2)2 - 4fl2 
The electric field emitted by the NV and cavity is described by 

^(+) = ^acNV + \fnaec + c.c, 



(F6) 



(F7) 



(F8) 



where e^v and Cc describe the spatial profiles of the emission from the NV center and the cavity. The emission is 
collected through a lens into an optical fiber (either single mode or multi-mode) with a set of orthonormal modes 
{/(^, cj)}. The field at the output of the fiber, which is directed at the detector, is given by 

E'^ = UriVl^eNv + V^ae^), (F9) 

where Up = ^% ^{'f{k,oj)){f{k,uj)-). The spectrum is obtained from the Quantum Regression theorem, S{uj) cx 
/-oo /o {E'^(t)E'~{t')) dtdt' . We make the assumption that the linewidth of the NV center is much broader than the 
NV-cavity detuning and all other loss or coupling rates: uj = ujc^ f^ ^ Id-td ^ Id-, ^- After normalizing the spectrum 
by the bare NV emission spectrum, we obtain 

S'{u) cc eNvUpeNV + 2-^[eNvUFece''^'^^/FW)--—^ r— ] + 

1 + l\UJ — ijJc)/ f^ 

eeUFeen^\-—T-^ r^P, (FIO) 

where F = {g{f^ fl))'^ /tzj denotes the Purcell factor at position r and orientation fl/\fi\ of the NV in the field of the 
cavity. In this experiment, it is not possible to calculate the term cnyUfCc directly because of uncertainty in the 
coupling of the NV and cavity modes to the fiber. In general, we therefore represent the detected spectrum as 

5rfH = C,+2C2W^^^W)-—-^ ^] + C3r(f)|-— — i ^\\ 

where the real coefficients Ci are obtained by fits to the data. We obtain good fits for the single mode fiber with 
C2/C1 ~ 0.6, and for the multi-mode fiber with C2/C1 ~ 0. 

Appendix G: Imaging of general samples 

The signal PL{uj^r) in Fig. 3 in the main text represents the convolution of the emitter's spatial profile, e(r), and 
the spectrometer signal given in Eq.l. Thus, PL{uj^ ^) — J ^{^~^)^d{^-, r')dl' ^ where dV is the displacement along the 
path of f . Since the NV center is less than 1 nm in size, we consider it as a delta function, so that PL{uj, r) = Sdioo^ r), 
giving the response function of our scanning probe. Once this response function is known, the SON maps a general 
sample e'{r) by first measuring the convolved image PL'{r) and then deconvolving by the known response Sd{oo^f') 
to obtain e'if). 
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